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In this paper we present a two-dimensional numerical simulation of a reconnection current layer 
in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large 
Lundquist numbers. We use realistic boundary conditions derived consistently from the outside 
magnetic field, and we also take into account the effect of the backpressure from flow into the 
the separatrix region. We find that within a few Alfven times the system reaches a steady state 
consistent with the Sweet-Parker model, even if the initial state is Petschek-like. 

PACS Numbers: 52.30. Jb, 96.60.Rd, 47.15. Cb. 
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Magnetic reconnection is of great interest in many 
space and laboratory plasmas , and has been studied 
extensively for more than four decades. The most im- 
portant question is that of the reconnection rate. The 
process of magnetic reconnection, is so complex, how- 
ever, that this question is still not completely resolved, 
even within the simplest possible canonical model: two- 
dimensional (2D) incompressible resistive magnetohydro- 
dynamics (MHD) with uniform resistivity r\ in the limit of 
S — > oo (where S = VaL/tj is the global Lundquist num- 
ber, L being the half- length of the reconnection layer). 
Historically, there were two drastically different estimates 
for the reconnection rate: the Sweet-Parker model 
gave a rather slow reconnection rate (Egp ~ S 1-1 ' 
while the Petschek |5j model gave any reconnection rate 
in the range from E$p up to the fast maximum Petschek 
rate -Epetschek ~ 1/ log'S 1 . Up until the present it was still 
unclear whether Petschek-like reconnection faster than 
Sweet-Parker reconnection is possible. Biskamp's simu- 
lations [|ll| are very persuasive that, in resistive MHD, 
the rate is generally that of Sweet-Parker. Still, his sim- 
ulations are for S in the range of a few thousand, and 
his boundary conditions are somewhat tailored to the re- 
connection rate he desires, the strength of the field and 
the length of layer adjusting to yield the Sweet-Parker 
rate. Thus, a more systematic boundary layer analysis is 
desirable to really settle the question. 

We believe that the methods developed in the present 
paper are rather universal and can be applied to a very 
broad class of reconnecting systems. However, for defi- 
niteness and clarity we keep in mind a particular global 
geometry presented in Fig. [I] (although we do not use 
it explicitly in our present analysis). This Figure shows 
the situation somewhere in the middle of the process of 
merging of two plasma cylinders. Regions I and H are 
ideal MHD regions: regions 1 represent unreconnected 
flux, and region II represents reconnected flux. The two 



regions I are separated by the very narrow reconnection 
current layer. Plasma from regions I enters the recon- 
nection layer and gets accelerated along the layer, finally 
entering the separatrix region between regions I and II. In 
general, both the reconnection layer and the separatrix 
region require resistive treatment. 
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FIG. 1. The global geometry. 

In the limit S — ► oo the reconnection rate is slow com- 
pared with the Alfven time ta = L/Va- Then one can 
break the whole problem into the global problem and 
the local problem || . The solution of the global problem 
is represented by a sequence of magnetostatic equilibria, 
while the solution of the local problem (concerning the 
narrow resistive reconnection layer and the separatrix re- 
gion) determines the reconnection rate. The role of the 
global problem is to give the general geometry of the re- 
connecting system, the position and the length of the re- 
connection layer and of the separatrix, and the boundary 
conditions for the local problem. These boundary con- 
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ditions are expressed in terms of the outside magnetic 
field B Vy0 (y), where y is the direction along the layer. 
In particular, B Vt o(y) provides the characteristic global 
scales: the half-length of the layer L, defined as the point 
where B y $(y) has minimum, and the global Alfven speed, 
defined as Va = B y ^(0)/y/Wp. 

In the present paper we study the local problem using 
the boundary conditions provided by our previous anal- 
ysis of the global problem (t| . Our main goal is to de- 
termine the internal structure of a steady state reconnec- 
tion current layer (i.e., to find the 2D profiles of plasma 
velocity and magnetic field), and the reconnection rate 
represented by the (uniform) electric field E. We as- 
sume incompressible resistive MHD with uniform resis- 
tivity. Perfect mirror symmetry is assumed with respect 
to both the x and y axes (see Fig. ^). 

This physical model is described by the following 
three steady state fluid equations: the incompressibil- 
ity condition, V • v = 0, the z component of Ohm's 
law, r)j z = E + [v x B] z , and the equation of motion, 
v • Vv = — Vp + [j z z x B] (with the density set to one). 

Now we take the crucial step in our analysis. We 
note that the reconnection problem is fundamentally a 
boundary layer problem, with S* -1 being the small pa- 
rameter. This allows us to perform a rescaling proce- 
dure [S inside the reconnection layer, to make rescaled 
resistivity equal to unity. We rescale distances and 
fields in the y-direction by the corresponding global val- 
ues (L, -Bo,j,(0), and Va), while rescaling distances and 
fields in the x-direction by the corresponding local val- 
ues: x -> x5 , v x -» v x V A S /L, B x -> B x By :0 (0)6 /L, 
E -> EB yfi (0)V A 5 a /L. Here, S = LS~ 1/2 is the Sweet- 
Parker thickness of the current layer. Thus, one can see 
that the small scale So emerges naturally. Then, using 
the small parameter Sq/L = S~ x l 2 <C 1, one obtains a 
simplified set of fluid equations for the rescaled dimcn- 
sionless quantities: 
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(where the first term on the right hand side (RHS) is the 
resistive term) and 



v • \7v y 



dy dx 



(3) 



In the last equation (representing the equation of mo- 
tion in the y-direction, along the current layer) the pres- 
sure term can be expressed in terms of B y (x,y) and the 
outside field Bo yV (y) by using the vertical pressure bal- 
ance (representing the ^-component of the equation of 
motion, across the current layer): 



p(x,y) = 



B io(y) 



B 2 y (x,y) 



(4) 



We believe that this rescaling procedure captures all 
the important dynamical features of the reconnection 
process. 

The problem is essentially two-dimensional, and re- 
quires a numerical approach. Therefore, we developed a 
numerical code for the main reconnection layer, supple- 
mented by another code for the separatrix region. The 
solution in the separatrix region is needed to provide the 
downstream boundary conditions for the main layer (see 
below). 

The steady state was achieved by following the true 
time evolution of the system starting with initial condi- 
tions discussed below. The time evolution was governed 
by two dynamical equations: 
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(Small artificial resistivity r\ y and viscosity v y were 
added for numerical stability.) The natural unit of time 
is the Alfven time ta — L/V A - The magnetic flux func- 
tion ^ is related to B via B x — —^f y , and B y = ^> x . 
At each time step, v x was obtained by integrating the 
incompressibility condition: v x (x,y) = — J Q X (dv y /dy)dx. 
Note that this means that we do not prescribe the incom- 
ing velocity, and hence the reconnection rate: the system 
itself determines how fast it wants to reconnect. 

We used the finite difference method with centered 
derivatives in x and y (second order accuracy). The time 
derivatives were one-sided. The numerical scheme was 
explicit in the y direction. In the x direction the resis- 
tive term d 2 ^/dx 2 was treated implicitly, while all other 
terms were treated explicitly. Calculations were carried 
out on a rectangular uniform grid. We considered only 
one quadrant because of symmetry (see Fig. ||). More 
details can be found in Ref. |9| . 
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FIG. 2. Computational box used in the numerical 
simulation. 
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The boundary conditions on the lower and left bound- 
aries were those of symmetry (see Fig. |^). On the 
upper (inflow) boundary x = xn m the boundary con- 
ditions were dv v /dx = (which worked better than 
v y = 0) and B y (x\i m ,y) = Bo, y (y) — the prescribed 
outside magnetic field. In our simulations we chose 
B Q , v (y) = B + (I - S ) v /l-y 2 with B = 0.3, con- 
sistent with the global analysis of our previous paper . 

The boundary conditions on the right (downstream) 
boundary cannot be given in a simple closed form. In- 
stead, they require matching with the solution in the 
separatrix region, which itself is just as complicated as 
the main layer. Therefore, we have developed a sup- 
plemental numerical procedure for the separatrix region. 
Noticing that in the separatrix region the resistive term 
should not qualitatively change the solution, we adopt 
a simplified ideal-MHD model for the separatrix. This 
model is expected to give a qualitatively correct picture 
of the dynamical influence of the separatrix region on 
the main layer, and thus a sufficiently reasonable down- 
stream boundary conditions for the main layer. In par- 
ticular, our model includes the effects of the backpressure 
that the separatrix exerts on the main layer. 

The advantages of our approach are: (i) use of the 
rescaled equations takes us directly into the realm of 
S — > oo; (ii) we do not prescribe the incoming veloc- 
ity v x (xn in , y) as a boundary condition: v x is determined 
not by the x-component of the equation of motion, but 
rather by v y via the incompressibility condition. As a re- 
sult, we do not prescribe the reconncction rate; (Hi) the 
use of true time evolution guarantees that the achieved 
steady state is two-dimensionally stable; (iv) we have a 
realistic variation of the outside magnetic field along the 
layer, with the endpoint L of the layer clearly defined as 
the point where Bo, y (y) has minimum (see Ref. 0). 

Let us now discuss the results of our simulations. We 
find that, after a period of a few Alfven times, the sys- 
tem reaches a Sweet-Parker-like steady state, indepen- 
dent of the initial configuration. In particular, when we 
start with a Petschek-like initial conditions (see Fig. |3^l), 
the high velocity flow rapidly sweeps away the trans- 
verse magnetic field B x (see Fig. ||). This is important, 
because, for a Petschek-like configuration to exist, the 
transverse component of the magnetic field on the mid- 
plane, B x (0,y), must be large enough to be able to sus- 
tain the Petschek shocks in the field reversal region. For 
this to happen, B x (0,y) has to rise rapidly with y in- 
side a very short diffusion region, y < y* -C L (in the 
case E init = 2Esp, presented in Fig. ||a, J/* = L/4), to 
reach a certain large value (B x = 2 for -Emit = 2i?sp) for 
y* < y < L. While the transverse magnetic flux is being 
swept away by the plasma flow, it is being regenerated by 
the merging of the B y field, but only at a certain rate and 
only on a global scale in the y-direction, related to the 
nonuniformity of the outside magnetic field B y fi(y), as 
discussed by Kulsrud |jj . As a result, the initial Petschek- 
like structure is destroyed, and the inflow of the magnetic 



flux through the upper boundary drops in a fraction of 
one Alfven time. Then, after a transient period, the sys- 
tem reaches a steady state consistent with the Sweet- 
Parker model. 
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FIG. 3. The current density j(x,y): (a) at t = 
for Petschek-like initial conditions with E^t = 2i?sp, 
and (b) in the final steady state (at t = 5), which 
corresponds to the Sweet-Parker solution. (All four 
quadrants are shown for clarity). 
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FIG. 4. The time evolution of the variation of the 
transverse magnetic field B x (0, y) along the midplane 
x — for the Petschek-like initial conditions. 

We believe that the fact that we rescaled x using the 
Sweet-Parker scaling does not mean that we prescribe 
the Sweet-Parker reconnection rate. Indeed, if the re- 
connecting system wanted to evolve towards Petschek's 
fast reconnection, it would then try to develop some new 
characteristic structures, e.g. Petschek-like shocks, which 
we would be able to see. Note that, if Petschek is cor- 
rect, then there should be a range of reconnection rates 
including those equal to any finite factor greater than 
one times the Sweet-Parker rate Egp- However, in our 
simulations we have demonstrated that there is only one 
stable solution and that it corresponds to E — Esp. In 
this sense we have demonstrated that Petschek must be 
wrong since reconnection can not even go a factor of two 
faster than Sweet-Parker, let alone almost the entire fac- 
tor of y/S. There seems no alternative to the conclusion 
that fast reconnection is impossible. 

It is interesting that in Petschek's original paper the 
length of the central diffusion region is an unde- 
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termined parameter, and the reconnection velocity v rcc 
depends on this parameter as VA{L/y*) 2 /VS. If y„ 
is taken as small as possible then Petschek finds that 
Vice ~ Va/ logCS 1 )" However, should be determined in- 
stead by balancing the generation of the transverse field 
B x against its loss by the Alfvenic flow (it should be re- 
marked that Petschek did not discuss the origin of this 
transverse field in his paper). As we discussed above, 
this balance yields y* ~ L, with the resulting unique rate 
equal to that of Sweet-Parker. This results are borne out 
by our time dependent numerical simulations. 

The final steady state solution is represented in Fig. ||b. 
It corresponds to a;u m = 5.0, yu m = 1.0, r\ y = v y — 0.01. 
We see that the solution is consistent with the Sweet- 
Parker picture of reconnection layer: the plasma param- 
eters change on the scale of order So in the x direction and 
on a global scale L in the y-direction. The reconnection 
rate in the steady state is surprisingly close to the typical 
Sweet-Parker reconnection rate E$p = ^^Vj^By o(0)- 
The solution is numerically robust: it does not depend 
on X\\ m , ?/ii m or on the small artificial resistivity r\ y and 
viscosity v y . 

Several things should be noted about this solution. 
First, j(x,y) — > (and B y {x,y) — > Bo iV (y)) monoton- 
ically as x — > oo, meaning that there is no flux pile-up. 
Second, as can be seen from Fig. Q B x (x = 0, y) ~ y 
near y = 0, contrary to the cubic behavior predicted by 
Priest-Cowley |l(J. This is due to the viscous boundary 
layer near the midplane x = and the resulting nonana- 
lytic behavior in the limit of zero viscosity, as explained 
m Ref. @. Third, there is a sharp change in B x and j 
near the downstream boundary y = y\ lm = 1, due to the 
fact that in the separatrix region we neglect the resistive 
term (which is in fact finite). 

It appears that the destruction of the initially-set- 
up Petschek-like configuration and its conversion into a 
Sweet-Parker layer happens so fast that it is determined 
by the dynamics in the main layer itself and by its inter- 
action with the upstream boundary conditions (scale of 
nonuniformity of Bq^), as outlined above. Therefore, the 
fact that our model for the separatrix does not describe 
flow in the separatrix completely accurately seems to be 
unimportant. However, for the solution of the problem 
to be really complete, a better job has to be done in de- 
scribing the separatrix dynamics, and, particularly, the 
dynamics in the very near vicinity of the endpoint of the 
reconnection layer. A proper consideration of the end- 
point can not be done in our rescaled variables, and a 
further rescaling of variables and matching is needed. 

To summarize, in this paper we present a definite so- 
lution to a particular clear-cut, mathematically consis- 
tent problem concerning the internal structure of the 
reconnection layer within the canonical framework (in- 
compressible 2D MHD with uniform resistivity) with the 



outside field Bo^ y (y) varying on the global scale along the 
layer. Petschek-like solutions are found to be unstable, 
and the system quickly evolves from them to the unique 
stable solution corresponding to the Sweet-Parker layer. 
The reconnection rate is equal to the (rather slow) Sweet- 
Parker reconnection rate, Egp ~ \/\/~S. This main result 
is consistent with the results of simulations by Biskamp 
JTT[ and also with the experimental results in the MRX 
experiment [ ^2j . 

Finally, because the Sweet-Parker model with classi- 
cal (Spitzer) resistivity is too slow to explain solar flares, 
one has to add new physics to the model, e.g., locally 
enhanced anomalous resistivity This should change the 
situation dramatically, and may even create a situation 
where a Petschek-like structure with fast reconnection is 
possible (see, for example, Refs. p3| , pT| ,p]| ) . 
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